Transport through a ID Mott-Hubbard insulator of finite length. 
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Transport through a ID Mott-Hubbard insulator of a finite length L is studied beyond pertur- 
bative approach. At special value of the low energy constant of the interaction we have mapped 
the problem onto the exactly solvable models and found current vs. voltage V at high temperature 
T > max(m,TL) and at low energy T,V < Tl (Tl ~ v c /L; v c : charge velocity). The result shows 
that for the strong interaction creating a large Mott-Hubbard gap 2m S> Tl inside the wire, the 
transport is suppressed near half-filling everywhere inside the gap except for an exponentially small 
region of V,T < T L exp(-2m/T L ). 
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Motivated by current experiments |l|] of Tarucha's 
group, a few works recently appeared where an effect 
of Umklapp processes on quantum transport through a 
ID wire has been studied. Previously, it has been al- 
ready shown by experiment Q and in theory |^-||] that 
there occurs no renormalization of the conductance due 
to the electron-electron interaction. Nonetheless, this 
renormalization has been observed Im. A recent calcula- 
tion Q based on Kubo formalism [^studied the effect of 
the Umklapp backscattering and concluded that it affects 
the zero temperature conductance. This contradicts the 
perturbative results of Ref. [p||lC|], where the quantized 
conductance G = e 2 /(■kTi) has been obtained at T = 
(below we use the units, where e = K = 1). Another 
point of view suggests jiT]|l| that in a non-perturbative 
case of large Mott-Hubbard gap 2m 3> v c /L = Tl {L: 
length of the wire; v c : charge velocity) the conductance 
is strongly suppressed at low energy inside the gap simi- 
lar to the band gap case. 

In this paper, to clear up the difference between 
the Mott-Hubbard insulator and the band gap one, we 
map the problem at low energies and at high temper- 
ature onto the exactly solvable models making use of 
a free fermion value of the constant g of the forward 
scattering inside the wire. The results are shown in 
Figs. 1 and 2. At low energies when T, V <C Tl 
(T: temperature; V: voltage), we have found that a 
new energy scale T x oc Tj,ea;p[— — [i 2 /Tl] ap- 
pears in the system if the chemical potential \l is small 
enough: y/ m 2 — /Tl ^> 1. Below T x the conduc- 
tance is not suppressed contrary to the expectation of 
Starykh and Maslov Jl2| and the current increases lin- 
early. Above this energy the current saturates and the 
conductance goes down as T x /T reaching small values 
~ exp[— 1\Jm 2 — (i 2 /Tl] at T ~ Tl- At high tempera- 
ture T> Ti,m we confirmed the asymptotical behavior 
of conductance: G = (1 — cst^(l + cosh^)" 1 )/^ pre- 
dicted before for Mott-Hubbard insulator in perturbative 
regime of small gap |P,^fj[ where the constant depends on 
m/T L . 

A brief physical explanation to these results follows. At 



low energies T < Tl and (i <C m the charge field is quan- 
tized inside the wire at its values related to the degener- 
ate sin-Gordon vacua. Rare low energy excitations tunnel 
through the wire with the amplitude oc exp[— tti/Tl] as 
(anti)solitons switching the quantized value of the field. 
The whole process of tunneling, however, includes trans- 
formation of the reservoir electron into the sin-Gordon 
quasiparticles and back. This transformation results in 
a non-trivial scaling dimension of the tunneling operator 
equal to 1 /2 for the Mott-Hubbard insulator connected to 
the Fermi liquid reservoirs independently of any parame- 
ters. In the case of the band insulator, this dimension is 
marginal (= 1): the transformation is trivial and does not 
introduce additional energy dependence. The infrared 
relevantness of the tunneling with the 1/2 dimension 
brings out above resonance at zero energy. Meanwhile, 
the exponentially small tunneling amplitude specifies the 
narrow width of this resonance equal to the crossover 
energy. Increase of \fi\ favours tunneling of the quasi- 
particles of the same sort and ultimately produces their 
finite density in the wire. Then the interaction between 
these quasiparticles described with the two-particle S- 
matrices ]l3| dependent on g emerges. At low momenta 
the S'-matrix for the quasiparticles of the same sort is 
inevitably free fermion like, as at g = 1/2. It manifests 
in the renormalization group (RG) flow derived from the 
Bethe anzats solution for the massive phase |l4|] of the 
sin-Gordon model and in the exponent calculated for the 
Tomonaga Luttinger liquid (TLL) phase at low density 
p"5[ . Increase of T, on the other hand, is expected to 
entail, first, a thermally activated behavior of the con- 
ductance ex exp[— 2m/T] at Tl < T < m |l(| and then 
a power law dependence at m < T. Since the effective 
value of <7, in general, scales with energy, the 1/T depen- 
dence we found for g = 1/2 may vary at higher energies 
r»m depending on the high energy value of g. 

Transport through the finite length wire under a 
constant voltage V between the left and right leads 
could be described in the inhomogeneous Tomonaga- 
Luttingcr liquid model (TLL) with the Lagrangian 
|,0: Jdx{J2 b £ b (x,cb b ,d t q> b )+£ bs (x,Vt,cb c ,cb s )}. The 
bosonic fields (f> b (x,t), b = c,s relate to the deviations of 



1 



the charge and spin densities from their average values 
as following: pb(x,t) — (d x (f>b(x,t))/(^/2ir), respectively. 
The first part of the Lagrangian describes a free electron 
motion modified by the forward scattering interaction. 
The second part of the Lagrangian introduces backscat- 
tering inside the wire. Only its term corresponding to 



the Umklapp process of four Fermi momenta transfer is 
important near half-filling. This term does not involve 
the spin field. Therefore, our consideration will be re- 
stricted to the charge field only. For the clean wire this 
field is characterized by the Lagrangian: 



where <p(x) — 9(x)9(L — x) specifies a one channel 
wire of the length L adiabatically attached to the leads 
x > L,x < and vf(Ef) denotes the Fermi veloc- 
ity(energy) in the channel. The parameter 4fep varies 
the chemical potential fi = 2kpv c = E t hr of the wire 
from its zero value at half-filling. Outside the Hubbard 
gap this parameter coincides with the momentum trans- 
fered by the backscattering: four Fermi momenta mi- 
nus a vector of the reciprocal lattice, and relates the 
present results to the perturbative ones J^]. The con- 
stant of the forward scattering varies from g c {x) = g 
inside the wire (x £ [0, L]) to g c (x) = goo = 1 inside 
the leads, and the Umklapp scattering of the strength U 
is introduced inside the wire. The charge velocity v c (x) 
changes from vf outside the wire to a some constant 
v c inside it. In the absence of the Umklapp scattering, 
v c ~ vpjg and < g < 1 is determined by the for- 
ward scattering amplitude of the bare short range inter- 
action between electrons. Approaching the half-filling 
put the Umklapp scattering on. It entails an essen- 
tial renormalization of the low energy value of g, which 
flows to its free fermion value g = 1/2 in the massive 
phase (lj] (fi < to) where the coefficient of the cos- 
term scales to ~ to 2 and on approaching this phase 
|l5f \fi\ \ to. This value of g will be assumed below. 
The zero frequency current through the wire equals I = 
V/ir+ < if, sc >, where the backscattering current is 

i bsc = -2E F U/v F J^dxsm{Ak F x + 2Vt + \f2(f> c {x)). It 
will be shown later that 2i:EfU is a half gap m, opened 
by the backscattering ([!]) in the charge mode spectrum 
inside the wire. 

I. High temperatures T > TL,m - The average backscat- 
tering current < Ib sc >= J D<fil}, sc exp{i J dt J dx(L c + 
£bsc)} can be written as a formal infinite se- 
ries in U. Each term of it is an integral 



of product of the free bosonic correlators < 
exp{i^/2(j)(x , t)}exp{— i>/2<j)(y, 0)} >. Such a corre- 
lator approaches its uniform TLL expression when 
x,L — x,y,L — y ^> v c /T. Substitution of this form 
into above series allowed us || to find L-proportional 
part of the backscattering current neglecting the bound- 
ary contribution in the perturbative case. However, 
the problem is not perturbative, in general, due to a 
finite gap 2to creation. Therefore, application of the 
uniform correlator will give us a part of the backscat- 
tering current oc min(L, v c /m) with the relative error 
0(max(Ti,,m)/T), which is of the order of ratio of the 
border piece oc v c /T to the essential part of the "bulk" 
one. This relates to the high-temperature asymptotics 
of the whole current. 
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FIG. 1. Schematic linear bias conductance G vs. tempera- 
ture near half-filling /j<Tl: curve 1 for the weak interaction 
m <^Tl, curve 2 for the strong one. 



Calculation of above series with the uniform TLL cor- 
relator is equivalent to expanding the value g = 1/2 into 
the leads. Following Luther and Emery |Q we map this 
bosonic Lagrangian onto the free massive fermion one 
0,(L2| with the density of Lagrangian 



E 

=R,L= 



i)+(d t ± v c d x )i> a - mip(x)[e^ x+2Vt ^+(x,t)ij R (x,t) + h.c] 



(2) 



Here i/jr(l) i s right (left) chiral fermion field. 
The fermionized backscattering current Ibsc = 
2im Jq dx[exp{i(4kF x + 2Vt)}t(i^(x,t)tpii(x,t) — h.c] 
is the doubled backscattering current for the fermions 



|L7J under doubled voltage. To find its average we just 
need to know the fermionic reflection coefficient R as a 
function of dimensionless energy u) — e/m: 
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R(co) 



(w 2 - 1) +sin 2 ( v / w 2 " 



(3) 



where = m/T^ denotes the dimensionless traver- 
sal time. The analytical continuation is assumed for 
\u)\ < 1. Since the chemical potential for the right/left 
chiral fermions is /i ± V, respectively, the total current 
can be expressed as 



V msinh(V/T) 



du> 



R(w) 



cosh((mw - fi)/T) + cosh(]//T) 



(4) 



Only the leading term in max(m, Tl)/T of the right hand 
side of (Q) is meaningful. Extracting it, we find the high- 
temperature asymptotics as following 



MV/T)B(mt L ) 



V m 

7 ~ 7cosh((/i)/T) + cosh(F/T) 
sm 2 {yfuj 2 ' 



(5) 



B(x) EE / duj 



Ix) 



(w 2 -l)+sin 2 ( v / w 2 



la;) 



where function B(x) increases as irx at small x > 
and approaches the constant ~ n at x 3> 
1. Accuracy of this calculation of (Q) may be 
written as a factor 1 + 0(max(m^Ti J )/max(T^V)) 
to (|) if < max(T, V) or as 1 + 



0([r7iaa;(r7i,T L )/maa;(r,l/)](maa;(r,\/)/^) 2 e^l/^), oth- 
erwise. The high-temperature conductance (Fig.|l|) 
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(6) 



approaches zero at T rs m if the gap is large enough 
tti/Tl ^> 1 and < m. 

2. Low energies T, V << Tl - To find a low energy model 
for our problem we have to integrate out all high energy 
modes. We will try to escape direct integration follow- 
ing Wiegmann's effective way of constructing the Bethe- 
ansatz solvable models [[l9| for the Kondo problem and 
for the screening of a resonant level |lj],|2(]]. First, let us 
substitute 4>/V2 instead of in (|l]). It makes fermions 
interacting inside the leads and non-interacting inside the 
wire. Their passage through the wire may be described 
with the one-electron S-matrix dependent on the elec- 
tron momentum. The interaction between electrons in 
the leads with some two-particle S'-matrix. Then the so- 
lution could be constructed if the proper commutation 
relations between the S-matrices are met. Being inter- 
ested in the variation of energy less then Tl around the 
Fermi level, one can simplify the solution keeping the 
one-particle S'-matrix constant equal to its value on the 
Fermi level. It leads us to the problem of one impurity 
in the TLL. 

For the weak backscattering, the Lagrangian of this 
problem can be written as 



dxCt 



2 X v 2 F 



d t 4>c{t, x) 



■in 



d x (j) c {t,x) 



4tt 



YT L u 
1 

TTVp 



i{2Vt + V2cj> c (t,Q)) 



(7) 



where we rescaled <fi back and introduced a new energy 
cut-off parameter YTl with dimensionless constant Y 
which will be specified later. Parameter u is related to 
the weak reflection coefficient as:w 2 = v F R(fi/m). For 
the strong backscattering the tunneling Hamiltonian ap- 



proach may be applied |2lJ] . It was associated [£5| to the 
dual representation using the field 9 mutually conjugated 
to cj> : [9 a (x), 4>(j{y)\ = i2irsgn(x — y). The appropriate 
Lagrangian reads 



dxCt 



2 V 



d t c (t,x) 



4-7T 



d x 6 c (t, x) 



4-7T 



YT L u' 

TTVp 



cos(Vi + 6 c (t, 0)/\/2) 



(8) 



with u' 2 = vp(l — R(n/m)) proportional to the free mas- 
sive fermion transmittance and the voltage multiplied by 
g factor Both these Lagrangian are, indeed, equiva- 
lent |2^1 if interaction dependent relation between u and 
u' is met The above model (0) or(§) character- 

izes the point scatterer of any backscattering strength at 
low energy ^5j. Although, the exact relation between u 
or v! and the bare parameters of the scatterer remains 
unknown. Our problem is dually symmetrical to that 
of Kane and Fisher: suppression of the direct current 
in their problem equals suppression of the backscatter- 
ing one in our case. This correspondence allows us to 



re- write their solution 23 2fJ as follows: 



7T 2 7T1 

T 

I = arctan(y/T x ), T = 
G = 



7T 

T, 



1 T 
ir z l 2 ttT 



(9) 



where ip denotes the digamma function and satisfies: 
ifj'(l/2) — 7r 2 /2, ip'(x) on 1/x, x — > oo, and a new en- 
ergy scale T x Hi varies from T x — YTlV^R at the weak 
backscattering^) to T x = YTl(1 — R)/tt at the strong 
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one (ra) . tive one ||[l0| . The latter was derived making use of the 

Let us, first, compare this result with the perturba- long-time asymptotics for the correlator: 

< e *.C*,t) e -*„(v.o) >= cst (aT L ) 29 ( . [ ^' TL l r ~— w d F{x)F{y) (10) 

v ' v sinh(TT(t — ia)l )sinh(Tr(—t + ia)l ) 



where a — 1/Ep and F was simplified: F(x) = 

(x/L) 9r U±,m=i( m±x / L ) sr2m±1 « caf x^L-asJ/L 2 ]^. 
r = (1— g)/(l+g) is the reflection coefficient of the charge 
wave scattering on the contacts of the wire with the leads. 
One can see that the current calculated with this form of 
the correlator will be consistent with the Lagrangian (Q) 
if we choose Y matching 

Y _ est (2n/T L y-^ 
sin(^/T L ) y/*r(l+r)J r+1/2 (n/T L ) [ ' 

This specifies Y as a constant of the order of 1 for 
\li\/Ti, < 1. Then the solution (||) coincides with the 
perturbative result |p|,|Io|, that is I — V/n oc — V 3 and 
G-1/ttoc-T 2 . 




FIG. 2. Schematic current I vs. voltage V near half-filling 
(i <C Tl'- curve 1 is zero temperature dependence, curve 2 is 
the high temperature T ^> Tl one, the dashed lines are the 
low voltage I = V/n and high voltage asymptotics. 



T for T < T x and approaches its exponentially small 
asymptotics G = T X /(4T) oc exp(-2m/T L )T L /T above 
T x < T < T L . As | /it | increases, the reflection coefficient 
i?(/i/m) on the Fermi level goes down and T x exceeds 
T L . finally, a pproachin g its weak backscattering value 
2T X = YT L /y/R(fi/m), where the perturbative consider- 
ation is applicable. 

In summary, we studied transport through a ID Mott- 
Hubbard insulator beyond perturbative approach. As- 
suming that 5 = 1/2 near the half- filling in agreement 
with the Bethe anztas solutions we mapped the prob- 
lem onto the exactly solvable models and found cur- 
rent vs. voltage at high temperature T > max(m,T L ) 
and at low energy T,V < Tl. The solution of these 
models shows, in particularly, that the high-temperature 
transport through the Mott-Hubbard insulator is simi- 
lar to the one through the band gap insulator at g — 
1/2. At low energies, however, there is always a regime 
where the transport remains unsuppressed in the ab- 
sence of the impurity backscattering. For the strong 
interaction resulting in the opening of the large Mott- 
Hubbard gap, the transport through the wire is sup- 
pressed near the half-filling almost everywhere inside the 
gap except for an exponentially small low energy region 
V,T <T L exp(~m/T L ). 
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Turning to the case tti/Tl > 1 we cannot use the per- 
turbative expression ([ll]) anymore: the perturbative se- 
ries is not convergent due to a finite gap creation. Then 
above non-perturbative consideration is necessary. Ap- 
plication of the solution (^|) in this case reveals a quite 
remarkable property of low energy transport through the 
Mott-Hubbard insulator. There is an exponentially small 
value of T x = (1 - R((i/m))YT L /n oc T L exp(-2m/T L ) 
for fi <C m. Hence, the zero temperature current / (Fig.|^) 
is not suppressed for the voltage less than T x and sat- 
urates at T x /2 value when T x < V < T L . Similarly, 
the conductance (Fig.|l|) displays a small decrease oc T 2 
below its zero temperature value l/7r with increase of 
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